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Abstract 

In this paper we provide some simple characterizations for the spher- 
ical harmonics coefficients of an isotropic random field on the sphere. 
The main result is a characterization of isotropic gaussian fields through 
independence of the coefficients of their development in spherical har- 
monics. 
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1 Introduction 



In recent years the study of real random fields on the sphere has received much 
attention, as this topic is a necessary tool in the statistical study of the CMB 
(Cosmic Microwave Background) radiation. The existing physical literature is 
huge, with particular emphasis on testing for Gaussianity and isotropy (e.g., 
for a small sample, [I], [TT]) as both these issues have deep implications 
for cosmological physics. See for reviews on this subject and where 
many more references can be found. We also mention j^J, pQ and [TU] for a 
mathematical treatment of the subject. 

A natural tool for this kind of enquiry is the development of the random 
field in a series of spherical harmonics. This raises some simple, but not so 
obvious questions. For an isotropic random field T, given its development in 
spherical harmonics (see §2), 



what properties must be expected to be satisfied by the coefficient a^ m ? Or, 
from a different point of view, which conditions must be verified by the coeffi- 
cients in order that the development above defines an isotropic random field? 

These questions are of interest both from a statistical point of view (for in- 
stance in order to devise a test suitable to detect non-Gaussianity or anisotropy) 
or from a probabilistic point of view (how to sample an isotropic random field). 

In this paper we provide some results in this perspective. More precisely, 
first we prove that, for an isotropic random field, the coefficients are necessarily 
uncorrelated. This fact is well known, but we give a proof that holds without 
the assumption that the field is mean square continuous. 

Then it is proved that each of the complex r.v.'s a^ m has a distribution 
whose phase is uniform on [— tt, tt]. This implies in particular that, for any 
isotropic field, the ratio Reae m / Ima,£ m is necessarily distributed accordingly 
to a Cauchy distribution. 

Finally, having remarked that if the field T is gaussian then the a^ m 's, i = 
0, 1, . . . , m — 0, . . . , I are independent, we prove that also the converse is true. 
Thus the only isotropic fields such that the % m 's, I — 0, 1, . . . , m — 0, . . . , £ are 
independent are those that are gaussian. This, which is the main result of this 
paper, is not a consequence of the central limit theorem, but results from some 





1=0 m=-l 
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classical characterization of gaussian random variables, through independence 
of some linear statistics. 

This result gives a rigorous proof of claims that can be occasionally found 
in the cosmological literature (see [3] e.g.) 

2 Isotropic random fields 

In this section we recall some well known facts about isotropic random fields 
T defined on the unit sphere S 2 = {x G R 3 : \x\ = 1}. We assume that 
these fields are isotropic in the strong sense, that is, their probability law is 
invariant with respect to the action of 50(3). The isotropy assumption can 
then be stated as follows: for all g G 50(3) and X\,...,x p G 5 2 , the two 
vectors 



have the same distribution. Throughout this paper by "isotropic" we mean 
isotropic in this sense. 

We shall use the spherical coordinates on S 2 i.e. x = ($, tp), where < $ < 
tt , < ip < 2tt. Also we shall use for g G 50(3) the parameterization through 
the Euler angles < a, 7 < 2tt and < (3 < n. Assuming the right-hand side 
exists, for each i = 1, 2, ... we can define the random vector 



Js 2 

where dx = sin $ d(p dd denotes the Lebesgue measure on S 2 and Ye, the vector 
of spherical harmonics defined by 



(T(gxi),...,T(gx p )) and (T(xi), . . . , T(x p )) , 




(2) 



Yt.(e,<p) 



Y tm (8,(p) 




P £m (cos^)e imv , for m > , 



Yi m (9, ip) 



(-l) m Y^ m (9,<f) ,form<0; 



(3) 



here Pi m (cos9) denotes the associated Legendre functions i.e. 



PUx) = {-iT{i-x 2 T' 2 —p t {x) , p t {x) 

m = 0,1,2, ...,£ , 1= 1,2,3, 



1 d e 




m\ dx* 
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A detailed discussion of the properties of the spherical harmonics can be found 
in Varshalovich, Moskalev and Khersonskii [T2] . chapter 5. Our purpose in this 
paper is to provide some characterizations of the probability law of the vector 
at under the isotropy assumption. It is a standard result that the functions 
Yt m m = —£, . . . , £ form a basis for the vector space of functions on S 2 which 
are restrictions of homogeneous harmonic polynomials of degree £. This vector 
space being invariant by the action of 5*0(3), for any £ and g G SO (3) there 
exist a (2£ + 1) x (2£ + 1) matrix D l (g) such that 

Ye.(gx) = D e (g)Y £ (x) , (4) 

This provides a representation of 50(3) on C 2e+1 which moreover is irreducible 
(see Vilenkin and Klymik |13j . §9.2.6 e.g.) The matrices D e (g) are the so-called 
Wigner's D-matrices whose entries are, in terms of the Euler angles 

D Lm> ( a > A 7) = e ~ ima dLm' (P) e" im ' 7 , 

where 

= (-iy- m 'i(£ + m)\(£-m)\(£ + m>)\(£-m')\] 1/2 



x 



mw(<-m,i-mO j3_y k -m-m' ^ gin §_yi-2k-m- 



k=0 



k\(£-m - k)\ {£ - m' - k)\ (m + m 1 + k)\ 



Note that d l mrn , (0) = 5™' , where 5™' denotes the Kronecker delta function. It 
is immediate that, the Lebesgue measure of S 2 being invariant by the action 
of 50(3), 



a? = I T(x)Yg (x) dx = I T(gx)Ye(x)dx 

Js* Js* (5) 

= / T(x)D\g- l )Y L (x)dx = D\g- l )a L , 
Js 2 

In particular, in coordinates, 

i 

aim = a tm'D e rn , m (g) . (6) 

m'=-£ 

It is well known that for an isotropic field which is continuous in mean square 
the development (JTJ) holds, the convergence being in L 2 (see jH] e.g.). 

It is also useful to point out that, because of (jHJ), the coefficients satisfy the 
following identity 

0>lm = ( — l) ma £,-m- 

In particular aeo is real. 
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3 General properties of the spherical harmon- 
ics coefficients 

It is well-known that, for mean square continuous and isotropic random fields, 
the spherical harmonics coefficients are orthogonal, i.e. 

%m,5g = < J ^C A , (7) 

the sequence {Ce}^ = \^,... denoting the angular power spectrum of these fields. 
For reasons of completeness we give now a proof of this fact. Actually our 
statement is slightly stronger. As usual we denote by A* the complex conjugate 
of the matrix A. 

Proposition 1 Assume T isotropic. Then 

a) for all £ such that E[\ae.\ 2 } < oo, 

Eai,a* L = C^e+i , 

where he+i denotes the (2£ + 1) x (2£ + 1) identity matrix 

b) forall£i,£ 2 such that -Z?[|a^.| 2 ] < oo, -E[|a^.| 2 ] < oo 

Ea £l a* e ^ = 

(in the sense of the {2£\ + 1) x (2£ 2 + 1) zero matrix). 

Proof a) Let us denote by the covariance matrix of the random vector at- 
Since the vectors and D (g)a^, have the same distribution, they have the 
same covariance matrix. This gives 

= D\g)V t D\gY = D\g)V,D\g)- 1 

Since D e is an irreducible representation of 50(3), by Schur lemma T e is of 
the form Qi^+i- 

b) The representations D £l and D l2 are not equivalent for £ x ^ £ 2 , having 
different dimensions. Therefore again by Schur lemma, the identity 

Ea ei at 2 , = D^{g)Ea h alD^{g)- 1 

can hold only if the right hand side is the zero matrix. 



5 



Remark 2 We stress that (for strongly isotropic fields) Proposition 1 is strictly 
stronger than the standard result on mean square random fields. Indeed, it 
is immediate to show that ET 2 < oo implies YaIi E\ai\ 2 < oo, on the other 
hand, it is not difficult to find examples where mean square continuity fails 
but the assumptions of Proposition f are fulfilled. Consider for instance the 
field: 



where bi 3rri3 = nai 3m3 ; we assume that the a£ imi 's (i = 1, 2, 3) satisfy Q whereas 
i] is a random variable with infinite variance (for instance a Cauchy). It is not 
difficult to see that the field T is properly defined and strictly isotropic; al- 
though (|SJ) is clearly an artificial model, some closely related field may be of 
interest for practical applications: for instance in CMB data analysis it is often 
the case that the observed field is a superposition of signal plus foreground con- 
tamination, and the latter may be characterized by heavy tails at the highest 
multipoles (point sources). In such cases, it is of an obvious statistical interest 
to know that the standard properties of the spherical harmonics coefficients 
still hold at least for the multipoles where foreground contamination is absent. 
It is immediate to see that ET 2 = oo, whence the field cannot be mean-square 
continuous; however (j2J is still properly defined for I = li, l 2 (simply exchange 
the integral with the finite sum), and therefore Proposition 1 holds for these 
two vectors of spherical harmonics coefficients. 

Our next proposition provides three further characterizations for the spherical 
harmonics coefficients of isotropic fields. 

Proposition 3 Let T be an isotropic random field (not necessarily mean 
square continuous); then for all a/, such that E\ai \ 2 < oo, 
a) for all m = 1, . . . , 2, 



b) for all m — 1, 2, ...,£, Re a^ m and Im a^ m are uncorrelated, with variance 
E(Rea £m ) 2 = E(lma £m ) 2 = C e /2. 

c) The marginal distribution of Rea^ m , Imo< m is always symmetric, that 
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Re aim = Im ae m and 




Cauchy 



is. 



Re a em = - Re a im , Im a tm = - Im a fc 
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Proof a) For (3 = 7 = 0, JBJ) becomes 

a £m — e a^m 

for all m = —£, . . . , £, < a < 2ir. This entails 

/ Rea£ m \ d_ I cos<^ siny? \ / Rea^ 
I Im a,£ m y I — sin ip cos </? y I Im a^ 

for all m = —£, ...,£, < ip < 2tt. Thus the vector *(Rea£ m , Imaf m ) has a 
distribution that is invariant by rotations, that is, in polar coordinates, they 
can be written in of the form 




J Rcos(6) (10) 

where R is a random variable with values in M + , whereas G is uniform in 
[— 7r, 7r]. This entails immediately that arctan(Rea^ m / Ima< m ) ~ Z7(— f , §); 
the result follows immediately. 

b) This property is well-known if T is mean square continuous. From (jlUj) . 

E[Re aim • im a^ m ] = / r 2 d[iR{r) I cos sin ^c/^ = . 

c) It suffices to take (po = ir in Q. 

Remark 4 It is interesting to note how Proposition ^ implies that no infor- 
mation can be derived on the statistical distribution of an isotropic random 
field by the marginal distribution function of the ratios (Rea£ m / Ima^ m ). On 
the other hand, it may be possible to use these ratios to implement statistical 
tests of the assumption of isotropy, an issue which has gained a remarkable 
empirical relevance after the first release of the WMAP data in February 2003 

It is clear that if the field T is gaussian, then the r.v.'s (ae m )e,m is a gaussian 
family. We prove now an independence result for this family of r.v.'s. Thanks 
to Prop osit ion [H these r.v.'s are uncorrelated, but one must be careful, since, in 
the case of complex r.v.'s, absence of correlation and joint gaussian distribution 
does not imply independence. 

Proposition 5 For an isotropic gaussian random field the r.v.'s a^ m , £ — 
0, 1, . . . , m — 0, . . . ,£ are independent. 



7 



Proof Let be (£,m) ^ (£',m f ), m > 0, m' > 0. Then a^ m is uncorrelated with 
both at' m i and <v _ m / = a^ m '. Thus 

^Kmttfm'] = 0, £/[a^ m a£' m /] = E[ae m ae> - m i] = 

and the statement follows from Lemma |H1 If one at least among m and m' is 
equal to 0, then the r.v. ae m (or ae> m > is real and independence follows from 
absence of correlation as for the real case. 

Lemma 6 Let Z\, Z 2 be complex r.v.'s, centered and jointly gaussian. Then 
they are independent if and only if 

E[Z 1 Z 2 ] = J E[Z X Z 2 ] = U (11) 

Proof In one direction the statement is obvious. Let us assume that (fTTj) are 
satisfied. Then, if we set Z^ = Xk + iYk, k = 1,2, then 

E[X X X 2 + Y 1 Y 2 ] + iE[-X 1 Y 2 + YiX 2 ] = 
E[X X X 2 - Y{Y 2 ] + %E\X X Y 2 + Y X X 2 \ = 

From these one obtains E[X X X 2 ] = 0, E[Y X Y 2 ] = 0, E[X X Y 2 } = and E\Y X X 2 ] = 
0. This means that each of the r.v.'s Xi,Yi, X 2 ,Y 2 is uncorrelated with the 
other ones, so that, being jointly gaussian, they are independent. 

Which is less obvious, is that the converse also holds. The following is the 
main result of this paper. 

Theorem 7 For an isotropic random field, let I be such that E\ai] 2 < 00. 
Then the coefficients (a^ , a-ei ■ ■ ■ , a>u) are independent if and only if they are 
gaussian. 

Proof We just need to prove the "only if" part. Fix mi > 0, m 2 > 0, so that 
the two complex r.v.'s a^ mi and a^ m2 are independent. Note that we are not 
assuming the independence of Re a^ mi and Im a^ mi or of Re a^ m2 and Im a^ m2 . 
Thanks to JSJ), the two vectors Q(. and D e (g)a^ have the same distribution. 
Thus the two r.v.'s 

1 t 
L i = X! De m', mi (9)ai m ' and L 2 = ^ D^, rri2 (g)a M 

m'=—£ m'=—l 
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having the same joint distribution as a £mi and a£ m2 , are independent. Fix g 
so that the angles a such that ma ^ kir for all integers m, k and j3 such that 
d mm (0) and d mni2 {(3) are different from zero for all m = —£, . . .,£ (note that 
such a (3 certainly exists, because the functions d^ (/3), m,m' = —£,...,£ 
are analytic and can vanish only at a finite number of values (3 G [0,71*]). For 
such a choice of g, thanks to ©, 



ReLi \ _ , e / cosmia sinmia \ / Rea^ 

ImLi I ^-^ mm i v I — sinmia cosmia / \ Ima^ 



and 



ReL 2 \ _ y\ ,t ,n\ ( cosm 2 a shir^a \ / Rea&, 
ImLx / ^ mm 2^ \ — sinm2a cos m,2a / V Ima^ 



where the 2x2 matrices on the right hand-sides are always full rank. By 
the Skitovich-Darmois theorem below (see Kagan, Rao and Linnik [7] e.g.), 
it follows that each of the vectors (Re a,£ m , Ima^ m ) is bivariate Gaussian; as 
(Re ai m , Im ai m ) are uncorrelated and have the same variance by Proposition 
[TJ then ai m = Re cti m + i Im ae m is complex Gaussian. 

Theorem 8 ( Skitovich-Darmois ) Let X\ , . . . , X r be mutually independent ran- 
dom vectors in R n . If the linear statistics 

r r 

are independent, for some real nonsingular nxn matrices Aj,Bj, j = 1, . . . , r, 
then each of the vectors Xx, . . . , X r is normally distributed. 

In particular Theorem[7|implies that if an isotropic random field is mean square 
continuous and the coefficients a^ m , £ — 0, 1, . . . , m — 0, . . .1 are independent, 
then it is gaussian. 

Remark 9 Proposition [7| shows that it is not possible to generate isotropic 
random fields by sampling non-Gaussian, independent complex-valued random 
variables ae m , m = —£, ...,£. This fact shows that, apart from the gaussian 
case, it is not easy to sample a random field by simulating the values of the 
random coefficients a^ m . In particular sampling independent values of the a^s 
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with distributions other than the gaussian gives not rise to an isotropic random 
field. 

We wish also to point out that it is indeed possible to construct a non- 
gaussian random field by choosing the random coefficients a,£ m , £ = 0, 1, . . . 
m = 0, . . . , I independent and with an arbitrary distribution. If they satisfy 
the conditions of Proposition Q and the series (JTJ) converges, they certainly 
define a random field on S 2 . But Theorem [7| states that such a field cannot 
be isotropic. In particular Theorem [7| does not follow by any means from the 
central limit theorem. 

Acknowledgment. The authors wish to thank Professor R. Varadarajan for 
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